Evolution of the nuclear modification factors with rapidity and centrality in d+Au 

collisions at ^/snn = 200 GeV 
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We report on a study of the transverse momentum dependence of nuclear modification factors 
RdAu for charged hadrons produced in deuteron + gold collisions at ^Jsnn = 200 GeV, as a function 
of collision centrality and of the pseudorapidity (jj = 0, 1, 2.2, 3.2) of the produced hadrons. We find a 
significant and systematic decrease of RdAu with increasing rapidity. The mid-rapidity enhancement 
and the forward rapidity suppression are more pronounced in central collisions relative to peripheral 
collisions. These results are relevant to the study of the possible onset of gluon saturation at energies 
reached at BNL RHIC. 

PACS numbers: 25.75.Dw, 13.18.Hd, 25.75.-q 



Studies of deep inelastic scattering of leptons on pro- 
tons and nuclei have revealed a large component of glu- 
ons with small-x (i.e. fraction of the nucleon momen- 
tum) that appears to diverge with decreasing x 0] . How- 
ever, it has also been suggested that the density of glu- 
ons remains finite due to the increased role of gluon- 
gluon correlations ('gluon fusion'), forcing an upper limit 
on the total number of highly delocalized small-x glu- 
ons 0, 0]. Phenomenological descriptions of DESY ep 
collider HERA and Fermilab results based on gluon 
saturation appear to successfully describe the data. Con- 
sequently, nuclei at high energies may be thought of as 
highly correlated systems of small-a; gluons. A QCD 
based theory for dense small-a; systems, termed the Color 
Glass Condensate (CGC) has been developed [||. 

Collisions between hadronic systems at a center-of- 
mass energy ^/sjvat = 200 GeV at the BNL Relativis- 
tic Heavy Ion Collider (RHIC) provide a window on the 
small-a; gluon distributions of swiftly moving nuclei. In 
particular, collisions between deuterons and gold nuclei 



in which hadrons with pt > lGeV/ c, mostly produced 
by quark-gluon interactions, are detected close to the 
deuteron beam direction, allow for probing the small-a; 
components of the wave function of the gold nuclei. It has 
been predicted that gluon saturation effects will manifest 
themselves as a suppression in the transverse momentum 
distribution below a value that sets the scale of the effect 
0, . The transverse momentum scale for the onset of 
gluon saturation depends on the gluon density and thus 
on the number of nucleons, and is connected with the 
rapidity y of measured particles by Q 2 S ~ A 1 / 3 e Ay , where 
A — 0.2 - 0.3 is obtained from fits to HERA data. Thus 
saturation effects are most evident at large y or pseu- 
dorapidity r), i.e. at small angles relative to the beam 
direction. At RHIC energies and at mid-rapidity the sat- 
uration scale for Au ions is expected to be ~ 2 GeV 2 

m 

We report measurements of transverse momentum 
spectra of hadrons from p+p and d+Au collisions at 
•J snn — 200 GeV in four pseudorapidity ranges around 
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FIG. 1: Bottom row: Invariant yield distributions for charged hadrons produced in d+Au and p+p collisions at y/s = 200 GeV 
at pseudorapidities 77 = 0, 1.0, 2.2, 3.2. Horizontal lines indicate bin width. Top row: Outline of the data samples in r\ vs pr 
collected with the BRAHMS spectrometers at the various angle and magnetic field settings (shifted in rf), gray lines are used 
to indicate the overlap of the settings. The dotted lines indicate the cuts applied in 77. 



77 = 0, 1, 2.2, and 3.2 (rapidities of the deuteron and 
gold nuclei are +5.4 and —5.4, respectively). A pion with 
transverse momentum p T = 2 GeV/ c, at these rapidities, 
probes the gluon distribution in the gold nuclei down to 
x values that range from 0.01 at 77 = to 4 x 10~ 4 at 
7/ = 3.2. We compare the yields from rf+Au collisions 
to those from p+p, scaled by the average number of bi- 
nary collisions (N co u) in a d+Au event. Results around 
mid-rapidity have previously been reported [3 E| • 

The data presented here were collected with the 
BRAHMS detector system ^(J , consisting of event char- 
acterization detectors and two rotatable magnetic spec- 
trometers: the Forward Spectrometer (FS) and the Mid- 
Rapidity Spectrometer (MRS). For the present studies 
the MRS was positioned at 90 and 40 degrees and the FS 
at 12 and 4 degrees with respect to the deuteron direc- 
tion. The experimental method and analysis techniques 
employed here are similar to those used for the study of 
d+ Au and Au+ Au collisions , except that the present 
data at 77 = 2.2 and 77 = 3.2 were analyzed using the front 
part of the FS detector systems only. The minimum bias 
trigger is estimated to select 91% ± 3% of the d+Au in- 
elastic cross section and 71% ± 5% of the total inelastic 
proton-proton cross section of 41 mb. The p+p yields 
have been corrected for trigger bias. Our trigger selects 
non-single-diffractive events and we estimate using the 
PYTHIA model that the correction should be 13 ± 5%, 
approximately independent of px and 77. 

Figure ^ shows the invariant yields of charged hadrons 
obtained from d+Au collisions and p+p collisions in nar- 
row pseudorapidity intervals around 77 = 0, 1, 2.2 and 3.2. 
The spectra at 77 = 2.2 and 3.2 are for negative hadrons 
only. Each distribution was constructed from indepen- 



dent measurements at several magnetic field settings, as 
shown in the upper panels of Fig. and is corrected 
for the spectrometer acceptance and tracking efficiency. 
The FS acceptance ranges from 2 to 4% and is known 
with an accuracy that ranges from 3 to 5% of those val- 
ues. No corrections for the finite momentum resolution, 
binning effects, absorption or weak decays have been ap- 
plied; the 15% systematic error on the spectra includes 
the contribution from these effects. However, at the over- 
lap between field settings at 4 degrees (between 1 and 2 
GeV/c), the systematic error is 20%. The momentum 
resolution of the spectrometers at the maximum mag- 
netic field setting is Sp/p = 0.0077p for the MRS and 
Sp/p = 0.0018p for the FS, (with p in GeV/c). 

Figure [21 compares the d+Au spectra to p+p distribu- 
tions using the nuclear modification factor defined by: 

1 d 2 N d+Au /dp T dr] 
KdAu = (N coll ) cPNZS/dprdr, ■ ( ) 

For the minimum-bias sample We estimate the mean 
number of binary collisions (N co u) = 7.2 ± 0.3, using 
the HIJING v. 1.383 event generator [TJ and a GEANT 
based Monte-Carlo simulation of the experiment. In 
these ratios most systematic errors cancel. Remaining 
systematic errors arising from variations in collision ver- 
tex distributions, trigger efficiencies and background con- 
ditions etc. are estimated to be less than 10% at 77 = 
and less than 15% at all other angle settings. 

Figure [2] reveals a clear variation of the RdAu as a 
function of pseudorapidity. At mid-rapidity, RdAu (v t > 
2GeV/c) > 1 shows a Cronin type enhancement |l2j 
as compared to the binary scaling limit. At 77 = 1 the 
Cronin peak is not present and at more forward rapidities 
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FIG. 2: Nuclear modification factor for charged hadrons at pseudorapidities r\ = 0,1.0,2.2,3.2. One standard deviation 
statistical errors are shown with error bars. Systematic errors are shown with shaded boxes with widths set by the bin sizes. 
The shaded band around unity indicates the estimated error on the normalization to {N co u}. Dashed lines at pr < 1.5 GeV/c 
show the normalized charged particle density ratio ^ JV 1 ^ aw/d^pp)^ • 
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FIG. 3: Central (full points) and semi-central (open points) R cp ratios (see text for details) at pseudorapidities r\ = 
0, 1.0, 2.2, 3.2. Systematic errors (~ 5%) are smaller than the symbols. 



(r/ = 3.2) the data show a suppression at all px- The val- 
ues of the RdAu ratios at low px are observed to be similar 
to the ratio of charged-particle pseudorapidity densities 
in rf+Au and p+p [l| collisions {n ] m) ^ffi^g 
shown in Fig. with dashed lines at px < 1.5GeV/ c . 

Figure |3 shows the ratio R cp of yields from collisions of 
a given centrality class (0-20% or 30-50%) to yields from 
more peripheral collisions (60-80%), scaled by the mean 
number of binary collisions in each sample. The central- 
ity selection is based on charged particle multiplicity in 
the range —2.2 < rj < 2.2 as described in |l3j. Since the 
peripheral collisions are similar to p+p, the R cp is domi- 
nated by the nuclear effects in the more central collisions, 
making the nuclear modification independent of the p+p 
reference spectrum. The data from the different central- 
ity classes are obtained from the same collider run. The 
ratios shown in Fig. are therefore largely free of sys- 
tematic errors associated with run-by-run Collider and 
detector performance, and wide i] bins can be used for 
each spectrometer setting. In contrast, the ratios shown 



in Fig. |21 must be constructed from two collider runs 
with different species. Smaller rj bins must then be used 
in order to include detailed acceptance corrections lead- 
ing to larger fluctuations. The dominant systematic error 
in the R cp ratios comes from the determination of (N co u) 
in the centrality bins. The shaded bands in Fig. [21 indi- 
cate the uncertainty in the calculation of (N co u) in the 
peripheral collisions (12%). We estimate the mean num- 
ber of binary collisions in the three centrality classes to 
be (iV c o 7 z 20% ) = 13.6 ± 0.3, (N™J 50% ) = 7.9 ± 0.4 and 

(K 



60-80%\ 
coll / 



= 3.3 ±0.4. 



There is a substantial change in R cp between rj = 
and the forward rapidities. At low pseudorapidity, the 
central-to-peripheral collisions ratio is larger than the 
semicentral-to-peripheral ratio, suggesting the increased 
role of Cronin like multiple scattering effects in the more 
violent collisions. Conversely, at forward pseudorapidi- 
ties the more central ratio is smallest indicating a sup- 
pression mechanism that depends on the centrality of 
the collision. In Fig. 01 we show R cp for the transverse 
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momentum interval pr = 2.5 — 4.0GeV/c. A fit to 
R C p ~ e ari yields a = -0.28 ± 0.03 for the central- 
to-peripheral ratio, a similar rj dependence as Q 2 S from 
HERA @, a = -0.13 ± 0.03 for the semicentral-to- 
peripheral ratio. 

The observed suppression of yield in d+Au collisions 
(as compared to p+p collisions) has been qualitatively 
predicted by several authors 0, 0, 0] within the 
framework of gluon saturation that includes the effects 
of 'quantum evolution' with rapidity. However, no de- 
tailed numerical predictions are yet available. These ap- 
proaches also predict the observed centrality dependence 
of the suppression at different pseudorapidities. Other 
authors pH . llflj have based their predictions on the two 
component microscopic HIJING model that includes a 
parametrization of perturbative QCD and string break- 
ing as a mechanism to account for soft coherent particle 
production, and 'gluon shadowing' as a method for reduc- 
ing the number of effective gluon gluon collisions. The 
HIJING model has been shown to give a good descrip- 
tion of the overall charged particle distribution in d+Au 
collisions 0, |i(J , and thus the \ow-px behavior of RdAu 
with nseudora.niditv. 
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FIG. 4: Evolution of the central/peripheral (full points) and 
semicentral/peripheral (open points) R cp ratios on pseudora- 
pidity. 

In summary, we observed a significant reduction of the 
yield of charged hadrons measured in c?+Au collisions, as 
compared to scaled p+p collisions at forward pseudora- 
pidities. This suppression for pr > 2GeV/c, is absent 
at mid-rapidity [jj EI , increases smoothly as the differ- 
ence in rapidity between the detected particles and the 
gold ion increases. Also, the change from mid- to for- 
ward rapidities is stronger for central collisions than for 
semicentral collisions, indicating a dependence on the ge- 
ometry of the collision. Such effects are consistent with 
the onset of saturation in the Au nuclei gluon density at 
small— a; values which modifies the shapes and magnitudes 



of RdAu and R cp at all transverse momenta. 

These results highlight opportunities for studying sat- 
uration phenomena in nuclei at RHIC. 
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